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Abstract 

In this work, we analyze the necessary conditions for a nonlinear AdS-type 
instability in the gravity-scalar field system. In particular, we discuss the nec¬ 
essary conditions for a cascade of energy to higher modes by applying results 
in KAM theory. Our analytical framework explains numerical observations 
of instability of spacetimes even when the spectrum is only asymptotically 
resonant, and the fact that a minimum field amplitude is needed to trigger 
it. Our framework can be applied for similar stability analyses of fields in 
(locally) asymptotically AdS spacetimes. We illustrate this with the example 
of the AdS soliton. Further, we conjecture on the possible reasons for quasi- 
periodic behaviour observed in perturbations of AdS spacetime. Certain 
properties of the eigenfunctions of the linear system dictate whether there 
will be localization in space leading to black hole formation. These properties 
are examined using the asymptotics of Jacobi polynomials. Finally we dis¬ 
cuss recent results on the AdS instability in Einstein-Gauss-Bonnet gravity 
in light of our work. 
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I Introduction 


Anti-de Sitter (AdS) spacetime is the arena for the AdS/CFT conjecture 
that provides a dual description of gravity in (an aymptotically) AdS space- 
time (with other helds) in terms of a non-gravitational held theory on its 
boundary. In this context, there is a huge body of work on held propagation 
in AdS spacetime treated in a linearized approximation. The importance of 
nonlinearities in AdS has since been realized due to recent work by Bizon 
and Rostworowski [1] that demonstrates that AdS spacetime is nonlinearly 
unstable to massless scalar held perturbations, however small. Their result 
is primarily numerical, and the end-point of the instability is the formation 
of a black hole in four dimensions and higher (see also |2])lfl They have 
demonstrated the nature of the instability by considering weakly nonlinear 
perturbations. The nonlinearity triggers a resonant instability wherein en¬ 
ergy is transferred to higher frequency modes (of the linearized system). A 
similar result holds for complex scalar helds [5]. 

In this paper, we investigate the necessary conditions for (i) an AdS'-like 
(resonant) instability to occur in held propagation in spacetimes, and (ii) 
given such an instability, the necessary condition for it to result in black hole 
formation. An aim is to pin down what features of held propagation in AdS 
such as a perfectly resonant spectrum of linear modes, or rehecting boundary 
conditions are crucial for the resonant instability. A more important reason 
is to investigate if asymptotically AdS spacetimes suher from nonlinear in¬ 
stabilities as AdS does. An obvious application of this is to the AdS/CFT 
correspondence. 

In section II, we show that a paradigm shift in viewing the AdS insta¬ 
bility is the key to understanding when resonant instability can occur in the 
gravity-held system. As we discuss in II.2, perturbation theory in AdS can 
be studied as perturbations of an appropriate integrable Hamiltonian — in- 
tegrability implies dynamics is restricted to be on tori in phase space, and 
instability due to the perturbation leads to destruction of these tori. The 
AdS instability is thus an instability in phase space. In II.3, we apply a 
theorem in nonlinear dynamics due to Benettin and Gallavotti [B] which is 
specihc to the integrable Hamiltonian that emerges in the gravity-scalar held 
system. Using this, we derive the necessary conditions for these tori to be 

^ In three dimensions, a black hole does not form, but nevertheless, initially smooth 
perturbations lose smoothness, and their higher Sobolev norms grow with time 0, i- 
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destroyed, and these in turn lead to conditions under which a instability in 
phase space can happen. If the operator governing linear perturbations is 
self-adjoint with the chosen boundary conditions, then the necessary condi¬ 
tion for instability is a linear spectrum that is discrete and either resonant or 
approximately resonant (in a precise sense). In the approximately resonant 
case, a minimum nonzero held amplitude (that depends on the deviation of 
the spectrum from perfect resonance) is needed to trigger the instability. An 
example is an asymptotically resonant spectrum (approaching resonant for 
high mode number) which arises in helds in a cavity in Minkowski spacetime 
with Neumann boundary conditions. Our analysis explains for the hrst time, 
the puzzling numerical studies which show in this case that the turbulent 
energy cascade can occur [7j, but only with a minimum nonzero scalar held 
amplitude |8] (see also [9]). In II.4 (and Appendix A), we demonstrate the 
application of our analysis to (locally) asymptotically AdS spacetimes by 
considering scalar held perturbations of the AdS soliton in weakly nonlinear 
perturbation theory. We identify a class of metric perturbations of the AdS 
soliton caused by the back-reaction of the held, for which the equations are 
very similar to the AdS scalar held system. Using our framework, we derive 
the necessary conditions for instabilities in this system. Further nonlinear 
aspects are discussed in section II.5. 

For the end-point of the instability to be a black hole, we also need a 
sufficient localization of energy in space. Given a resonant instability that 
transfers energy to higher modes, in section III, we argue that a necessary 
condition for localization of energy is localization of the high mode linear 
eigenfunctions, at least for small held amplitudes. This is demonstrated 
explicitly for AdS, where the asymptotics (high mode behaviour) of the mode 
functions (Jacobi polynomials) is studied using the Mehler-Heine formula 
for various dimensions. Localization of higher mode functions is least for 
AdSs consistent with the observation in [3] that black hole formation does 
not occur in AdS^ despite a very fast transfer of energy to higher modes 
(in comparison with higher dimensions). This leads to an understanding of 
when an energy cascade can lead to black hole formation in an asymptotically 
AdS spacetime. Section II and III together derive necessary conditions for an 
AdF-like nonlinear (in)stability in asymptotically AdS spacetimes — this can 
now be checked by a simple investigation of the spectrum and eigenfunctions 
of the linear problem. 

Bizon and Rostworwski’s work has led to several efforts to obtain a bet¬ 
ter analytical understanding of the instability. Among these has been the 
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observation m that the AdS instability is similar to the evolution of a sys¬ 
tem of coupled nonlinear oscillators hrst studied by Fermi, Pasta, Ulam and 
Tsingou (FPUT) [11]. Similar to the FPUT problem, the authors of [10] 
separate the fast and slow time scales in the system. A careful analysis of 
this separation using rigorous renormalization group methods was done in 
|12j . The two-time framework has been used since in [13], [E]. This has 
uncovered quasi-periodic and recurrent behaviour in the gravity-scalar held 
system for times 0(l/e^), e being the held amplitude at initial time. Our 
work, on the other hand, gives predictions about the behaviour of the system 
for exponentially long times in 1/e^. We also conjecture on reasons for this 
quasi-periodic behaviour in section II.5. Another area of investigation in this 
topic is whether special initial perturbations are stable. In particular, hne- 
tuned periodic perturbations which were constructed in [15] and boson stars 
seem to be stable (see also Ha)- These are the nonlinear versions of 

single modes of the linearized analysis (‘oscillons’). For evolution of initial 
perturbations comprising of two modes of the linearized analysis, numerical 
data is inconclusive, but this has led to lively debate [IB], [H], and a more 
recent numerical analysis is in [2D]- It has been surmised that there may 
be quasi-periodic solutions consisting of many modes for special initial con¬ 
ditions [20]. The issue of whether gravitational interaction localizes energy 
in position space has been investigated in [21] by using shells of helds. Our 
work adopts a different approach — of studying asymptotics of eigenfunc¬ 
tions which explains the dimension dependence of localization, such as no 
horizon formation in AdS^. 

The results of our paper are summarized in section IV. In section V, we 
have an extended discussion of open problems in this area. In particular, we 
discuss a recent result in Einstein-Gauss-Bonnet gravity coupled to a scalar 
held of amplitude e [22] which shows complex behaviour in comparison to 
the pure gravity case. We argue that under a held rescaling these equations 
are identical to the pure gravity case to O(e^), but not beyond. The details 
of this computation are given in Appendix B. 
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II Necessary conditions for AdS-type resonant 
instability in gravity-scalar field system 

In this section, we obtain a deeper understanding of the AdS instability using 
nonlinear dynamics. This paradigm shift enables us to understand precisely 
when held perturbations in a spacetime can lead to a resonant instability 
with a cascade of energy to higher frequencies as observed in AdS. We obtain 
necessary conditions for such a resonant instability in the gravity-scalar held 
system, and this analysis can be directly applied to scalar (and other) helds 
in asymptotically AdS spacetimes. 


II. 1 Summary of gravity-scalar field set-up in AdS 

We hrst summarize the Einstein gravity-scalar held set-up used by Bizon 
and Rostworowski to uncover the instability of AdS spacetime [T]. Their 
hrst paper considered (3 -|- 1) dimensions, but for future use, we write the 
equations in (d-|-l)-dimensions. Let A = The coupled held equations 

are: 

Gab + Agab = SttG (^daddbcj) - ^gabiddf^ . (II. 1) 


g^’^A/aA/bd = 0. (II.2) 

The scalar held and metric are assumed to have spherical symmetry with the 
metric ansatz being 


ds^ = 




COS^ X 


'—Ae "^^dP + A ^dx^ + sin^ xdSP 


d-l) 


(II.3) 


where is the standard metric on the {d — l)-sphere. The range of 

coordinates are —oo < t < oo and 0 < a; < n/2. A, 6 and 0 are functions of 
t and X only. The pure AdS solution is given by A = 1, 5 = 0 and 0 = 0. In 
this system, rehecting boundary conditions are imposed on the scalar held0 
A solution to the Einstein-scalar held system in weakly nonlinear per¬ 
turbation theory is obtained by the expansion in a small parameter e > 0 


"'^For technical subtleties associated to the equations of the gravity-scalar field system 
with this boundary condition, see |23) . 
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as follows (to make contact with previous literature, we use units: SnG = 
d — 1; c = 1): 

OO OD OD 

= 02,+16^^+'; A = 1 - 5^ 6 = ■ (11-4) 

j=0 j=l j=l 

The zeroth order (in e) background metric is AdS with 0 = 0. The hrst-order 
(in e) equation is 0, + L0i = 0 where 

L = -- — ][zi-dx [(tan^^-^a;)^,,] . (II.5) 

tclll X 

The solution obeying reflecting boundary conditions at the AdS boundary 

X = 7r/2 is of the form X] 0ii = X] cos(a;jt + /3j)ej{x). ej{x) are the 

j j 

eigenfunctions of the linear operator L with eigenvalues Uj = {2j + d)^. L is 
self-adjoint on the space of square integrable functions with inner product 

/•7r/2 

<f,g>= / f{x)g{x)taLn'^~^xdx. 

Jo 

The (orthogonal) eigenfunctions are 

i \ 2'\/j!(j -I- d — 1)! ^ ^tt 

ej[x)= -cos a:P/ ^ (cos 2a:); (II.6) 

where are the Jacobi polynomials and j are non-negative integers. 

At order e^, as a result of the back-reaction of the scalar held, we have: 

da = - / cos(?/)sin(|/)[(0;)2 + (0i)^]d|/. (II.7) 

do 

, sinxcosa: , 

Aa = -— —I — / [(0i) + (0i) ]tan'^ d (H.S) 

t cLil X J Q 

At the third order an inhomogeneous equation is obtained 

03 + -^03 = >S'(Aa, da, 0i) = ~20i(da + Aa) — 0i(da J- Aa) — 0i(A2 + d^). (11.9) 

We can expand 03 = as {ej{x)} form a basis on the space of 

j 

functions on which L is self-adjoint (assuming 03 lies in this function space). 
Projecting fill.91) onto the basis {ej{x)} results in 

Cj + u'^Cj = (5, ej{x)) (II. 10) 
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which is an infinite set of forced harmonic oscillator eqnations for the co¬ 
efficients Cj(t) = ( 03 ,ej(x)). Triads of freqnencies in {S,ej{x)) of the form 
Uj^} can lead to resonances when ±Uj = This is dne 

to the perfectly resonant linear spectrnm of the scalar field in AdS. Not all 
secnlar terms can be removed, leading to a resonant instability and transfer 
of energy to high freqnency modes over time. 


II.2 Instability of the scalar field-gravity system from 
the nonlinear dynamics point of view 


We now characterize the resonant instability in a different way, which yields 
many benefits. Consider the eqnation for the scalar field fill.2D . with the 
metric flll.Sp solving the Einstein eqnation flll.ip . The compntations in the 
previons section show that we can expand the scalar field as 

in terms of the basis {ej{x)} of eigenfnnctions of a self-adjoint operatorjfl 
Then, in weakly nonlinear pertnrbation theory flII.4D . dj{t) can be expanded 
as a series of odd powers of e. At order e, the dj{t) are jnst harmonic, with 
freqnencies uj. The next (third) order correction is given by solntions of 
a harmonic oscillator eqnation with forcing terms from the nonlinearities, 
as in flll.iop . The eqnations governing dj{t) thus arise from a Hamiltonian 
that can be viewed as a perturbation of linear harmonic oscillators, where 
the perturbation couples modes, leading to forcing terms at higher orders. 
In principle, the gravity-scalar field system has an infinite number of such 
coupled nonlinear oscillators. But at any finite time (however large), we may 
expect significant energy transfer to only a finite number of modes. In studies 
of the AdS instability, most numerical simulations show a decaying power 
law behaviour of energy per mode with mode frequency, close to black hole 
formation times [2]. In a two-time formalism, conservation laws have been 
found [1^ (which may hold at times of the order ^) that lead to inverse 
energy cascades [26]. Thus we may view the Hamiltonian for dj{t) as a 
perturbed system of n linear (undamped) harmonic oscillators, where n can 
be large but is finite. As we shall see, the analysis that follows from this 


® Existence results for the spherically symmetric Einstein-scalar field system [24], [25] 
suggest that in the case of even spacetime dimensions, the basis for the scalar field of 
the full nonlinear system may be a corrected version of the linear eigenfunctions with 
corrections being higher order in e. This will not change the discussion of this subsection 
and only add one more term to the source in (III.101) . 
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explains puzzling numerical results in the literature. 

The Hamiltonian for dj{t) can be written as an integrable part (n linear 
harmonic oscillators) plus a perturbation of order e^. By the Liouville-Arnold 
theorem, the integrable system of n linear oscillators can be described in 
terms of action-angle variables {(/*, 0*)} G xT”, i = and dynamics 

takes place on n-tori, which are level surfaces of the action variables. We 
use the notation (/, 6) to collectively describe the 2n variables, dropping 
subscripts. Denoting the Hamiltonian of the integrable n linear oscillators 
by ho{I), the Hamiltonian of the perturbed system is ho{I) + e‘^f{I,9). We 
will now analyze this system without making any specific assumptions about 
which in our system gives rise to source terms seen, for example, in 
flH.lOp . f{I,9) is in general not integrable, but we can apply procedures in 
Hamiltonian perturbation theory (HPT) to this near-integrable Hamiltonian, 
to hnd a canonical transformation that makes it approximately integrable to 
some order in e. This ensures dynamics still takes place on some possibly 
perturbed version of the n torus of the linear system —we characterize this 
as stability of the system. However, as is known, there is an obstruction 
to hnding such a canonical transformation due to the small denominator 
problem in HPT. This implies that if the frequencies of the linear oscillators 
are resonant, with ratios of frequencies rational, even a small perturbation 
can destroy the corresponding tori of the linear system and dynamics is not 
restricted any more to be on them. As the system samples more of the 
phase space, we expect energy to be transferred to a large number of modes, 
via resonances. In AdS, the frequencies of the linear modes are resonant. 
Therefore, this is another way of viewing the resonant instability of AdS, as 
an instability in phase space. 

In fact, frequencies which are near-resonant (ratios close to rational) can 
also lead to the small denominator problem and lack of integrability. If this 
occurs, a consequence will again be destruction of tori and a nonlinear trans¬ 
fer of energy across modes, leading to an instability. Such frequencies do 
occur in the gravity-scalar held system and instabilites have been observed 
numerically in such cases, as further evidence for our picture. One example is 
the analysis of a massless scalar held in a cavity in Minkowski spacetime with 
Neumann boundary conditions on the boundary of the cavity [7j. An insta¬ 
bility was found numerically in this case. The spectrum is not resonant, but 
for large frequencies, becomes asymptotically resonant (ratios of frequencies 
become close to rational for higher modes). Further, it has been observed 
that a certain minimum amplitude is required to trigger instability in this 




case, as opposed to arbitrarily small perturbations in the fully resonant case 
[H]. Massive fields in a cavity in Minkowski spacetime (with either Dirich- 
let or Neumann boundary conditions) have a nonresonant spectrum that is 
asymptotically resonant. Instability has been found numerically in this case 
for small held amplitudes by Okawa, Cardoso and Pani [9]. These results 
have been puzzling so far due to the view that instability is associated with 
a resonant spectrum. Viewing the instability as destruction of tori in phase 
space provides a reason for instability even when spectra are only asymptot¬ 
ically resonant. We will show further that it also explains the details and 
nuances of numerically observed features. 


II.3 Applying the Benettin-Gallavotti theorem to un¬ 
derstand numerical results 


We will now apply a theorem by Benettin and Gallavotti [ 6 ] in nonlinear 
dynamics (KAM theory), to the system of n perturbed linear harmonic oscil¬ 
lators with hamiltonian H{I,6) = ho(/) + This represents weakly 

nonlinear perturbations of the gravity-scalar held system. This theorem gives 
the dynamics of the perturbed system over long times which are an expo¬ 
nential of a power of far longer than the time scales ~ ^ studied in the 
two-time formalism. We would like to understand the reasons for features 
observed numerically in the case of asymptotically resonant spectra. 

Let Ui,i = 1,..., n denote the frequencies of the linear oscillators (modes of 
the scalar held in the linearized analysis) and let oj G denote the frequency 
vector with components Ui. Let k G 2’"' — {0} denote a vector of integers. 
The condition for a resonance is that a; ■ k = 0 for some k G — {0}. If 
the entire spectrum or part of it is perfectly resonant, we can hnd n linearly 
independent k satisfying the resonance condition in the fully resonant case, 
and rii < n in the case where a subspace of the frequency space of dimension 
rii is resonant. Let us now consider a spectrum that is not perfectly resonant 
(either fully or partly). This can be quantihed by a Diophantine condition 
on u for all k E — { 0 } and some 7 > 0 , namely 0 : 


\(jj ■ k 



( 11 . 11 ) 


\k\ denoting supremum (over i) of \ki\. Clearly by choosing large integers 
for ki we can get arbitrarily close to the resonance condition. 7 quantihes 

®For a discussion of the Diophantine condition, see p.217, m- 
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how close the frequencies are to being perfectly resonant — asymptotically 
resonant or nearly resonant spectra will satisfy this condition only for small 

7- 

Benettin-Gallavotti theorem: 

Consider H{I, 6) = ho(/) + 9). Let the frequencies obey the Diophan- 

tine condition fill. lip for some 7 . Further, denote by 1/F^ the timescale for 
oscillation of the linearized system, E = ||a;||. With some further analyticity 
conditions on Benettin and Gallavotti showed the following: 

For initial data (/(O), 0(0)), and for amplitude of the perturbation < Cq, 
where Cq = {^)'y and D > 0 is a constant determined by n and analyticity 
parameters of /, [B], [25] : 


||J(f)-J(0)||<C(-); 

<^0 


( 11 , 12 ) 


valid for times t such that 




-fiAh*’ 


(11.13) 


where C is a constant, b = ... . 

4(n+l) 

In the absence of a perturbation, the action variables are constant, and 
flII.12p indicates that in the presence of the perturbation, the action variables 
remain close to the initial values for exponentially long times in comparison 
to the time scale of oscillations of the linear system This is thus an 
exponential stability result in the case when the frequencies are nonresonant. 
It is valid for a perturbation of size where eo depends on 7 in fill.lip . 

The stability result will be of practical use if the frequencies are strongly 
nonresonant, i.e., 7 is not very small. It is important to note that when the 
condition < Cq is not fulhlled, it does not imply automatically that there 
is an instability. Rather, the above stability result does not hold and due to 
the small denominator problem, an instability is very likely. We will discuss 
the caveats to this statement in II.5. 

If the frequencies are asymptotically resonant (i.e., approaching a reso¬ 
nant spectrum as j —)■ cx)), 7 will be very small, and so will Cq. When the 
spectrum is asymptotically resonant, the instability, if present, will need a 
minimum nonzero amplitude of the order of Cq to trigger it, in contrast to 
what was observed in the fully resonant case. However, that minimum am¬ 
plitude required for triggering it can be very small if the initial condition 


10 










has higher frequency mode components for which 7 is small (leading to small 
eo)0 Also, the weight of the higher modes with frequency Uj in the initial 
condition matters. We must have the effective amplitude of the higher modes 
of the held greater than Cq for triggering instability. If the perturbation has 
this minimum amplitude, we expect the time scales for the instability in the 
asymptotically resonant case to be similar to that in the fully resonant case 
as both are a result of lack of (approximate) integrability of the perturbed 
system. This explains the numerical observations in 0. 0 and [9]. For a 
massless scalar held in a cavity of radius R in Minkowski spacetime ana¬ 
lyzed by Maliborski [7j, with Dirichlet boundary conditions on the cavity, 
the spectrum is resonant. However, with Neumann boundary conditions, an 
asymptotically resonant spectrum is obtained where Uj are the roots of the 
equation tan{ujR) = cUjR. For large j, Uj ~ |(j + |) - ^{j + + C»(^)- 

Even for j ~ 100, retaining only the hrst (resonant) piece is a good approx¬ 
imation. This implies that if the initial condition has modes of this order, 
by our discussion, the minimum amplitude required for instability will be 
small. Okawa, Cardoso and Pani observe [9] that even with Dirichlet bound¬ 
ary conditions in the cavity, there is a nonresonant spectrum if the scalar 

held is massive. The frequencies are given by Uj = \J where iR 
is the scalar held mass. This is an asymptotically resonant spectrum when 
<< 1- Then, for sufficiently large j satisfying this inequality, we will 

have Uj ^ ^ + + 0(^). Even if fiR = 10 (at the higher end of the range 

considered in 0 ). for j ~ 100, retaining the hrst (resonant) term is a good 
approximation. If the initial condition has modes of this order or higher, an 
instability can occur with small amplitudes of the held, as shown in [9]. The 
normalization used in [9] was G = 1, and diherent from that used in our 
paper. Keeping k = SttG explicitly in the equations only implies that the 
magnitude of the perturbation term in the Hamiltonian is | e^, and k = Sir 
in [9]. 

A review of analytic aspects of the AdS instability, including our work on 
applications of nonlinear dynamics results to asymptotically resonant spectra 
can be found in [5U] . 

A note of caution is in order here. These results have been arrived at for 
perturbations of a hnite number n of linear oscillators. Simply taking n ^ 00 

^It has also been argued in [55] that higher modes need to be excited for an instability 
in this case. 
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is physically meaningless, because all oscillators are not on an equal footing. 
In fact, numerics shows a decaying power law for energy per mode with mode 
frequency. There are inhnite dimensional versions of KAM theory, almost 
all of which have assumptions such that effectively only a hnite number of 
degrees of freedom are signihcant and the results are similar to what we have 
mentioned. 


II.4 Fields in (locally) asymptotically AdS spacetimes 
— the AdS soliton 

We now illustrate, by means of an example — the AdS soliton — the applica¬ 
tion of the framework of the previous subsections to many (locally) asymptot¬ 
ically AdS spacetimes. Let us consider, for simplicity, the gravity- massless 
scalar held system in weakly nonlinear perturbation theory where the scalar 
held and the metric are expanded as: 


9ij — 9{o)ij + 9(2)ij + •••• ; <h — e<hi -|- -|-. (11-14) 

We now consider the case when g{o)ij is a locally asymptotically AdS 
metric which is static, diagonal, and with no horizons or singularities. An 
example of this is the AdS soliton metric (written in hve dimensions): 


r 


ds^ = 


P 


1 - 


dA -|- (dx^)^ -|- {dxA — dt 


+ 1 - 


(11.15) 

The coordinate r is restricted to r > tq and r is a periodic coordinate 
which must have a period /5 = ^ to avoid a conical singularity at r = tq. 
For the detailed construction of this metric, we refer to [31] . 

The Klein-Gordon equation for $1 in a static, diagonal background fi'{o)p, 
reduces to the form 


and can be analyzed by separation of variables into modes. For the AdS 
soliton, let the scalar held depend only on r and t. Then <Fi = 0i(r) exp{iu}t) 
and L = —:^dr[{A — r^Adr]. The mode frequencies obey the eigenvalue 
equation L(f)i = This equation cannot be solved exactly, but the large 

mode frequencies have been obtained in various ways, either numerically or 
in a WKB approximation [32], [33], [33], [33], [33]- Here, we give the WKB 
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result [20] valid for /r^ > 5 (we have set the dependence of their ansatz on 

coordinates and to zero, and their calculation of glueball mass spectra 
then gives the squared frequency): 

a;^~?7.(n + l)—^—hO(n°); (11.16) 

where n are large positive integers. We observe that this is an asymptotically 
resonant spectrum as for large u, ca ~ Cn where C is a constant. The nor¬ 
malizable eigenfunctions corresponding to these frequencies can be evaluated 
asymptotically, and vanish at the boundary. They are also regular aX r = R. 

Now, for the gravity-scalar held system, with <f)(r, f), let us consider met¬ 
rics Qij in the class 

ds^ = 

P f{r) 

(11.17) 

where /(r) = r^(l — ((j). 

This class preserves the planar and circular symmetries of the background. 
It has been considered in the context of nonlinear gravitational perturbations 
of the soliton in [37]. In Appendix A, we give the full Einstein-scalar held 
equations for this system. The AdS soliton background corresponds to the 
functions a, b,c,d = 0. fill.141) gives the expansion for 4) as well as the metric 
functions as 

a = 026 ^ -|- a^e'^ . 

and similarly for b, c, d. Considering the equations in Appendix A at order 
yields coupled equations for the functions 02 , C 25 <^ 2 . Interestingly, at order 

we get 

—rPdtdt^s + (5r^ — rQ)dr^3 + P — = 

rP{-d2 + &2 + C 2 2d2)dt^i ri‘^{-2a2 + 2b2)dtdt^i 

_(^5 _ rgr)(a 2 — b '2 + Cg + 2d'2)dr^i (11.18) 

This equation can be cast in a form very similar to what was obtained in 
firm for AdS; 

dtdt^s + L4>3 = S{a2, &2, C 2 , ^ 2 , *^* 1 ) (11.19) 

where L is the operator governing the linear modes and S = —(—02 + 62 + 
C2 -|- 2 d 2 )dt^i -\- -— ^4 (^2 ~ ^2 A C2 2 d 2 )dj.^i 2(02 — b 2 )dtdt^i. 
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The simple structure of the equations order by order in e observed in 
the scalar held -AdS system is reproduced for this example as well, by con¬ 
sidering appropriate classes of perturbations. We can stretch the analogy 
further by considering the class of reduced perturbations a = b = 0. For 
this class, the second order equations can be decoupled to obtain C 2 and d 2 
explicitly analogous to flll.Sp in the AdS case. The details of this calcula¬ 
tion are given in Appendix A. If the eigenfunctions of the linear operator L 
provide a basis for expanding $ 3 , then by projecting (III.191) onto this basis, 
we obtain the perturbed linear harmonic oscillator structure discussed exten¬ 
sively in this section. Of course, a rigorous existence/uniqueness analysis for 
the Einstein-scalar held equations for this class of perturbations as done for 
AdS in [21] would be complicated in this case. Thus we do not know if the 
eigenfunctions of L are the appropriate basis, order by order, for <F. Nev¬ 
ertheless, this exercise shows that it is possible to apply the entire analysis 
of the previous subsections to static, diagonal, asymptotically AdS spaces, 
with suitably restricted perturbations. In particular, if the eigenfunctions of 
L are an appropriate basis, then the asymptotically resonant spectrum of the 
AdS soliton implies that a necessary condition for instability is presence of 
the higher modes in the initial condition and a perturbation of a minimum 
nonzero amplitude. As is well-known, the AdS soliton is the minimum of 
energy in its asymptotic class [3l]. The instability, if it exists, would be for 
(hnite) scalar held excitations of this background. 

We note that a recent paper considering scalar held perturbations of 
quasi-periodic solutions in AdS [38| has also obtained an asymptotically res¬ 
onant spectrum for the linear modes of the held. 

II.5 Discussion on further nonlinear aspects 

The KAM (Kolmogorov-Arnold-Moser) theorem [31] is often quoted in the 
AdS instability literature as relevant for long-time behaviour of the gravity- 
scalar held system. Unlike the Benettin-Gallavotti result (on stability for 
exponentially long times), with the same assumptions, the KAM theorem 
guarantees stability as t —>■ cx) provided the unperturbed Hamiltonian ho(/) 
obeys the non-degeneracy condition 7 ^ 0 . A set of linear harmonic 

oscillators do not obey this condition as this determinant is zero. Therefore, 
we do not see any straightforward application of this theorem, unless the 
perturbed Hamiltonian can be rewritten as a nonlinear Hamiltonian obeying 
this, plus a small perturbation. The Nekhoroshev theorem [10] which shows 
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exponential stability even in the case of a resonant spectrum, with some 
assumptions (a geometric condition on the Hamiltonian, and the nondegen¬ 
eracy condition) also does not seem to apply here. 

The striking similarity of the AdS gravity scalar held system to coupled 
nonlinear oscillator Hamiltonian problems such as the FPUT problem m has 
been used in HDl to uncover quasi-periodic behaviour at times of the order 
Quasiperiodic solutions in various contexts in the gravity-scalar held system 
have been found by [I5], [26], [16]. Selection rules that restrict energy transfer 
to higher modes (and ‘inverse energy cascades’) have been found in [13], [12], 
[26] and more recently, [TT]. Similar quasiperiodic behaviour/selection rules 
have been seen in a system of gravity and self-interacting scalar helds [TT] . 
[H], [13], [S]. At hrst glance, this is puzzling, since we generically expect a 
resonant instability in AdS, due to its spectrum. 

We conjecture that the reason for quasiperiodic behaviour at times of the 
order ^ could be the ‘closeness’ of the nonlinear scalar held-gravity system 
to a nontrivial integrable model (i.e., not just a set of linear oscillators) to 
some order in e. The FPUT system of coupled nonlinear oscillators can also 
be viewed as a perturbation of linear oscillators. However, this does not 
explain its complex behaviour. It has been shown that the FPUT model 
displays quasiperiodic behaviour for certain time scales which can be stud¬ 
ied in a two-time framework (such as that employed in some of the papers 
mentioned). At one scale, FPUT dynamics is indistinguishable from the in¬ 
tegrable Toda model m, m (see also [U] for the connection of the FPUT 
model to the KdV equation). In the Toda model, coefficients corresponding 
to nonlinear mode mixing are simply zero, thus those resonant instability 
channels are absent. At the second time scale, the FPUT model can be 
regarded as a perturbation of the integrable Toda model (rather than as a 
perturbation of linear harmonic oscillators). This accounts for the quasi¬ 
periodic behaviour for long times (after which one may expect a departure 
from Toda-like behaviour). The relation of the FPUT model to the KdV 
equation m also suggests that these quasiperiodic solutions are the ana¬ 
logues of the KdV breathers. We conjecture that a similar explanation holds 
for the quasi-per iodic behaviour of the gravity-scalar held system. There 
may be several ways of writing the perturbed H{I,6) = ho{I) + e^f{I,6). 
We have taken ho{I) to be the Hamiltonian for a system of linear harmonic 
oscillators, but if there is an integrable model with Hamiltonian ho{I) ‘close’ 
to our nonlinear system, H{I,9) = ho{I) + e^f{I,9) could be viewed as a 
perturbed integrable model. The identihcation of the integrable model ho{I) 
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is a hard problem, but one that we expect will open up exciting avenues of 
future study. 


Ill Asymptotics of Jacobi polynomials and 
significance for localization of field in space 

In the last section, we have seen some necessary conditions for resonant 
transfer of energy to higher frequencies. However, even if transfer of energy 
takes place to higher modes, we will still need sufficient localization in space 
for black hole formation. The key to understanding localization in the AdS- 
scalar held system is the asymptotics of the linear eigenfunctions for high 
frequency. We argue that given a resonant instability in the gravity-held 
system in a spacetime, localization in space of high frequency linear modes is 
a necessary condition for black hole formation, at least when held amplitudes 
are small. 

In studies of the AdS instability, most numerical simulations show a de¬ 
caying power law behaviour of energy per mode with mode frequency [2]. 
However, at sufficiently large times (close to black hole formation), it is rea¬ 
sonable to expect that there would be signihcant energy transfer to some 
high frequency modes. We will discuss the asymptotic form of the modes of 
the scalar held in AdS for high frequency in various dimensions. In fact, we 
have checked that the asymptotic (Mehler-Heine) formula we will discuss in 
this section for the value of the Jacobi polynomial at the maximum x = 0 
only dihers by about 1% from its actual value for mode number j = 100. 
In ^^54, even for j = 10 there is only an error of about 4%. The results of 
this section will be valid for mode numbers in this range as long as sufficient 
energy transfer has taken place to them. We expect their spatial behaviour 
to be useful in investigating the localization of the held in space. 

The normalized eigenfunctions of the scalar held in a hxed AdSd+i back¬ 
ground are given by 



(HI.20) 


where 



(HI.21) 
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Now, we want the asymptotic behaviour of ej{x) for large j. There are 
two regimes, one when x G ( 0 , 7 r/ 2 ) and the other at end-points x = 0 and 
X = ti/2 . First, let us consider x G ( 0 , 7 r/ 2 ). For large j, from the Darboux 
formula, we have | 18 ] 


(--1 

2 

j 


'(cos 2 a;) ~ cos{2jx + (j)) 


J'- 


( 111 . 22 ) 


Here, 0 = and {p)j = p{p+l)...{p+j — 1) refers to the Pochhammer 

symbol. 

A{x) = 2-^{smxp^-^'>/^{cosxp^+^^/\ 

Using this in (1III.20P and using Stirling’s formula, we obtain that for large j, 


ej{x) 


2 ^ [cos( 2 ja;) cos 0 — sin( 2 ja;) sin 0 ] 

\/F (sina;)U-i)/2(cosa;)U+i)/2 


(III.23) 


This formula is not valid at x = 0 and x = 7 r/ 2 , however we can certainly 
consider x —?• 7r/2, and ej{x) approaches zero in this limit for all d > 2. 
Similarly, let us consider the regime when j is large and x small, such that 
jx is large. In this regime, behaviour varies with dimension: 

(i) AdSs, or d = 2: 0 = —Clearly, for small x, ^ is the envelope function 
of rapidly oscillating ej{x). The envelope function increases slowly as x —)■ 0. 
For small x. 


ej{x) 


[Cl cos( 2 jx) — C 2 sin( 2 jx)] 


(ii) AdS^, or d = 3: 0 = —7r/2. For small x, ej{x) ~ We know 

that this function localizes near x = 0 and is in fact proportional to the delta 
function d(x) as j —)• 00 . 

(iii) For AdSd+i with d > 3, the enveloping function of ej{x) is as 

X —?• 0 and therefore localizes even better near x = 0. 

We now inspect the values of the normalized eigenfunctions at x = 0 
using the Mehler-Heine formula for Jacobi polynomials, which gives their 
values in a neighbourhood of x = 0 and x = ti/2 for large j |1H]. From the 
formula, we can easily check that ej{TT/2) = 0. Dehne 2 ; by 4 = 2x. Then, in 
a neighbourhood of x = 0, 

lim j-"P"’^(cos -) = (z/2)-"J„(z). (III.24) 

j^OO J 
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Ja{z) is the Bessel function. In the regime where j is large, x is small, such 
that z is large, the Mehler-Heine formula gives the same asymptotics as the 
Darboux formula flIIL23D . 

If i is large and a; —)■ 0, such that 2 ; —)■ 0, then Ja{z) 

Here, a = ^ — 1. This can also 


Thus, for X = 0, ~ 

be seen from the relation of the Jacobi polynomial to the hypergeometric 
function. We next want 6^(0). From flIII.20D . we need to estimate the large 

i behaviour of the normalization factor • Using Stirling’s formula 


for the Gamma function, this is y U Thus 6^(0) ~ r{a+i) ^ • This increases 

as j increases for all d > 2, with the increase being the smallest for d = 2, 
i.e., AdS^. 

We thus conclude that localization of the eigenfunctions around x = 0 for 
large j increases with dimension d. 

Thus if energy is transferred to higher frequency modes for which the 
asymptotic Darboux and Mehler-Heine formulae apply, their spatial behaviour 
and localization properties lead to localization of the held, at least for a wide 
variety of initial conditions with a single pulse such as the Gaussian in [T]. 
We see that localization increases with dimension, and is least for d = 2 
(AdSs). In AdSs, numerical studies show a very fast transfer of energy to 
higher frequencies (in comparison to higher dimensions), yet no black hole 
formation 0 . ffl- We speculate that at least part of the reason may be that 
the eigenfunctions localize the least in three dimensions. For asymptotically 
AdS spacetimes, a similar study of the eigenmodes at high frequency needs 
to be done to see if localization happens. Localization of eigenfunctions for 
high frequencies is a necessary (but not sufficient) condition for an instability 
to result in black hole formation for a wide variety of initial conditions. 

Guriously, the asymptotic Darboux formula for the scalar held mode in 
AdS^ is similar to the exact form of the mode eigenfunction for the scalar 
held in a cavity of radius R in 4-dimensional Minkowski spacetime with 
Dirichlet boundary conditions. Both when the scalar held is massless [7] and 
when it is massive [9], the eigenfunctions corresponding to the linear modes 

are ej{r) = _ Glearly for Uj large, these eigenfunctions localize 

(approaching a delta function in the limit Uj —)■ cxd). There is an instability 
in this system leading to horizon formation, consistent with our arguments. 

We would like to end the section with a word of caution. If one were 
to start with an initial condition with multiple pulses at diherent locations. 
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then the interactions between them could be very complicated [21] • In such 
a scenario, one cannot conclude naively about localization from the eigen¬ 
functions alone. One could also have initial data comprising of a few modes, 
but with an amplitude so large that black hole formation occurs without the 
need for resonant transfer of energy. Further, in gravity with a Gauss-Bonnet 
term [22], this term results (in some cases) in a pulse breaking up into mul¬ 
tiple pulses, with relative phases between themselves, and again asymptotics 
of eigenfunctions may not be useful. Asymptotics of eigenfunctions for x ^ 0 
have also been used recently in a different context, to compute interaction 
coefficients between modes in [19|, [50]. 

IV Summary 

The AdS instability makes it imperative to identify the necessary conditions 
for instabilities to be caused by nonlinearities in the gravity-field system. 
In particular, one would like to know when arbitrarily small field perturba¬ 
tions can cause an instability. This is difficult to probe numerically and an 
analytical understanding is needed. 

In this paper, we provide, for the hrst time, an analytical framework to 
identify the necessary conditions for an AdF-type instability in a gravity- 
scalar held system. It also predicts the late-time behaviour of arbitrarily 
small AdS perturbations. In this framework, the AdS instability is under¬ 
stood as a dynamical instability of the scalar held in phase space. Dynamics 
in the linearized approximation is integrable, and nonlinearities can lead to 
breakdown of integrability. If this happens, the scalar held can sample a 
higher-dimensional subset of phase space than the integrable case, with the 
only restriction coming from energy conservation. This manifests as nonlin¬ 
ear energy transfer across modes, leading to an instability. 

We apply a theorem in nonlinear dynamics, from KAM theory, to de¬ 
rive the necessary condition for such an instability. The spectrum of the 
linear modes of the held is crucial. If the spectrum is strongly nonreso¬ 
nant (i.e, ratios of frequencies are irrational numbers badly approximated 
by rational numbers in a precise sense), the gravity-scalar held system is 
stable for long times (even with rehecting boundary conditions). Thus, a 
necessary condition for instability is a spectrum that is resonant or approx¬ 
imately/asymptotically resonant. This is not sufficient, but an instability is 
likely in this case. The reason is a breakdown of integrability in the nonlin- 
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ear system due to the small denominator problem in nonlinear dynamics. A 
perfectly resonant spectrum, in general, leads to a resonant instability for ar¬ 
bitrarily small field amplitudes. If the spectrum is asymptotically resonant, 
we show that such modes can trigger an instability if they are excited by 
the initial condition, and if the perturbation has a certain minimum nonzero 
amplitude (whose value depends in a precise sense on the deviation of the 
spectrum from perfect resonance). This explains puzzling numerical results 
for a scalar held with an asymptotically resonant spectrum in a cavity in 
Minkowski spacetime 0 . 0 , 0 . 

Our framework is applicable to (scalar, gravitational and other) helds in 
many asymptotically AdS spacetimes for a wide class of non-dissipative (self- 
adjoint) boundary conditions. As an example, we consider the AdS soliton- 
scalar held system in hve dimensions, and identify a class of perturbations 
for which the structure of the equations bears a close resemblance to the AdS 
case. Since the spectrum of the linear modes of the held in the AdS soliton 
background is known to be asymptotically resonant, we predict the necessary 
condition for an instability leading to an energy cascade to higher modes. 

The outcome of a resonant instability in a gravitational system need not 
be a black hole. The necessary (but not sufficient) condition for that to 
happen given a resonant instability is a strong localization of modes of the 
linear system for higher frequencies. We have argued that in AdS^+i, this 
happens for a wide range of initial conditions due to the asymptotics and 
localization properties of Jacobi polynomials. Localization is least in AdS^. 
This provides the best explanation for lack of black hole formation in AdS^ 
despite efficient energy transfer to higher modes. 

These results together allow us to address the issue of stability of asymp¬ 
totically AdS spacetimes with respect to scalar helds and gravitational per¬ 
turbations (which have been discussed in [HI], [2S])- The next direction is 
a systematic study of asymptotically AdS spacetimes both with respect to 
their linear spectrum (under scalar and gravitational perturbations) and the 
localization properties of their eigenfunctions. This relatively simple exer¬ 
cise will tell us about their stability or the outcome of an instability. In the 
next section and in Appendix B, we have also offered a reason for surprising 
differences in the gravity-scalar held system between Einstein gravity and 
Einstein-Gauss-Bonnet gravity observed in |22] . 
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V Discussion, future directions 

As discussed in section II, fascinating quasi-periodic behaviour has been ob¬ 
served in numerical studies of the AdS instability problem until times t ~ ^ 
for some initial conditions. We have conjectured, based on our framework, 
that the reason could be proximity of the gravity-scalar held Hamiltonian to 
a nonlinear integrable Hamiltonian to 0{e^). This is similar to the Fermi- 
Pasta-Ulam system, whose quasi-periodic behaviour is explained by the fact 
that it is a perturbation of the integrable Toda model. A very interesting 
problem is the identihcation of the nonlinear integrable Hamiltonian that is 
close to the gravity-scalar held Hamiltonian to 0{e^). The various selection 
rules and results from the two-time formalism for studying the AdS instabil¬ 
ity will be crucial to solving this problem. 

One class of spacetimes to which our results do not apply are space- 
times with horizons, such as AdS black holes. This is due to the fact that 
quasinormal mode boundary conditions on perturbations at the horizon are 
dissipative, and the results we have used from nonlinear dynamics do not 
apply to systems with dissipation. The issue of nonlinear stability for some 
rotating black holes under gravitational perturbations has been discussed in 
|29]. However, it is known that even in the linearized approximation, small 
Kern-Ads' black holes are unstable due to superradiance isg. Numerical 
work on linear evolution of scalar helds in small AdS black hole spacetimes 
in the time domain conhrms this [5l]- It is very important to understand 
how the nonlinear effects impact the end-point of this instability. Recently, 
the set-up for full numerical nonlinear evolution of a scalar held in asymptot¬ 
ically AdS black hole spacetimes has been developed |S^. The results from 
this program may lead to insights on the role of nonlinearities in black hole 
perturbations. We are actively working, from an analytical point of view, on 
applications of nonlinear dynamics of dissipative systems to these exciting 
problems. 

A potentially rich and complex topic is the study of nonlinear instabilities 
in the gravity-scalar held system in theories of gravity with higher curvature 
corrections. The gravity-scalar held system was investigated recently in hve 
dimensional Gauss-Bonnet gravity [22] for the same set-up as that of Bi- 
zon and Rostworowski (reviewed in H.l). There are surprising diherences to 
Einstein gravity — in particular, the Gauss-Bonnet term leads to greater dis¬ 
persion of the scalar held. One observation in [22] is particularly interesting: 
for e (held amplitude) less than a certain nonzero critical value (correspond- 
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ing to the critical mass for no black hole to form), there is not even a resonant 
transfer of energy to higher scales. This needs to be analytically understood; 
however, we offer a reason for why the dynamics with the Gauss-Bonnet term 
is different from the pure Einstein case. Let us denote the Gauss-Bonnet cou¬ 
pling by A 3 . The equations of the Einstein-Gauss-Bonnet (EGB) scalar held 
system are given by flVIII.44p and flVIII.dSp in Appendix B. Studying the 
system order by order in weakly nonlinear perturbation theory, there are 
corrections to the corresponding equations in Einstein gravity beyond the 
linear order in e. We have found that, by dehning a rescaled held V’ by 

^ _ / 

(1 - 2A3)V2 

the equations for the EGB scalar held system for (j) are identical to the Ein¬ 
stein gravity-scalar held equations up to O(e^) for ip. The Gauss-Bonnet 
term cannot be scaled away at order and higher. The detailed equations 
up to O(e^) before and after rescaling can be found in Appendix B. Inter¬ 
estingly, precisely the same rescaling has been observed before in a two-time 
formalism (which is valid to 0(e^)) in |56j^ . The authors of [56] interpret this 
observation holographically using the AdS/GFT correspondence and conjec¬ 
ture that nonequilibrium dynamics of the conformal held theory (GET) dual 
in the presence of the Gauss-Bonnet term is similar to that of the GET in 
the absence of it up to times 0 (l/e^). 

Thus, at least at order e^, a perturbation in Gauss-Bonnet gravity is 
dynamically equivalent to a larger held perturbation in Einstein gravity — 
larger by the factor This factor can be signihcant if A 3 is large. 

Even if e were small, this factor being large would paradoxically imply that 
the system is dynamically equivalent to Einstein gravity to order but in 
a regime when weakly nonlinear perturbation theory is not expected to be 
valid. We therefore expect complex behaviour. As well, if e is not very small, 
clearly the perturbative picture breaks down and the Gauss-Bonnet term 
can play an important role as it introduces new nonlinearities. The regime 
that would be interesting to study is one where e and A 3 are both small. 
In this regime, at every order and higher, the terms involving the Gauss- 
Bonnet coupling A 3 are not the leading terms growing in time. So a priori, 
it does not seem that there would be drastic diherences to Einstein gravity 

®Our conventions for the cosmological constant differ from those in this paper, but after 
taking this into account, the rescaling factors are the same. 


22 










until late times when the perturbation has grown, and the Gauss-Bonnet 
term could dominate. The conjecture in [SS] is precisely the holographic 
dual of this observation. Indeed, as the authors of [56] point out, holography 
constrains the value of A 3 to be small, and this seems to be the physically 
interesting regime. This issue merits further numerical investigation, and 
a deeper analytical understanding. We hope to report on these exciting 
problems in the near future. 
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VII Appendix A 

For the class of hve dimensional metrics given by (III.1711 and a massless scalar 
held <F(r, f), the Einstein-scalar held equations can be written as {k = 8nG): 

Rab + 'j29ab = ^ab (VII.25) 


1 




0 


where Wg b = Kj Tgb - \gabT) and Tgb = d g^db^ - \ggb{,d^Y. 
Plugging flll.lTp into flVII.2511 & flVII.2611 we get 


(VII.26) 


.-26 r 




e^“r/'(l -|- ra') + e ^“/{2 -|- — rb' + rd + 2rd' + Ara' — r'^a'b' 


+ raid + 2r^a'd! + r^a!'} + e^^r{db + dc — b^ — c‘^ — 2d^ -|- 2dd — b 


— c — 2d} 




Wu 


(VII.27) 
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-.—2a r 


— — e^°‘rf'{?> + ra' — rb' + 3rc' + 2rd') — 2{e^“r/(ra'^ + 2a' 


2r^f 

— ra'b' + rc''^ + Ad' + 2rd'‘^ — 3b' — rb'c' — 2rb'd' + ra" + rc" + 2rd") + 

4e2b 


e^^l^{db — be — 2bd — 6 ^ — 6 )} 


+ 


/ 


= Wr^ 


(VII.28) 


^-2a-2b+2c 


f 


2 /V^ 


e^"r^f" + e^“r/'(3 + ra' — rb' + 3rc' + 2rd') + 2{e^“r/(rc' 


/2 


+ 3c' + ra'c' — rb'd + 2rc'd' + rc") + e^^d{dc — be — 2cd — — c)} 

Ae^^f 


+ 




= W^ 


(VII.29) 


— e 


— 2a—2b+2d 


e^“r/'(l + rd!) + e ^“/{2 + rd + Srrf' + r'^c'd' + 2r‘^d'‘^ + ra' 


+ r^a'd' - rb' - r%'d' + r^d"} + e^H^{dd -bd-cd- 2d^ - d) 


+ Ar'^e^'^ = (VII.30) 

with Wtt = K<i>2, Wrr = & Wrr = = 0 . 


— rde~°‘~^^ (dtdt^ + (—a + 6 + c + 2(i)5t$ j — ror)5r5r$+ 

(5r^ — rQ)c}r<h + (r® — rQr){a' — b' + d + 2d')drA>^ = 0 (VII.31) 

Let us now consider a reduced class of perturbations with a = b = 0. Putting 
a = 6 = 0 in the equations flVII.27tlVII.30p . one obtains: 


rf 

d 



rfd 

d 


+ 


2rfd' 

d 


-c^ -2d'^ -d-2d-A 




Wtt 


(VII.32) 


1 

2rf 


rf" + 3/' + Srf'd + 2f'rd' + 2f{rc'^ + Ad' + 2rd'^ + rc" + 2rd") 


+ 


4 

7 


= IV„ 


(VII.33) 
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-1 


r^f" + r/'(3 + 3rc' + 2rd') + 2r/(rc'^ + 3c' + 2rc'd' + rc") 
+ 4 = 0. 


2 r2 
+ 4c(i + 2c) 


l\2c^ 

(VII.34) 


— rf — r‘^f'd' — 2/ — bfrd' — 2fr^d'‘^ — frd — fr^c'd' — fr^d" + l*{cd 
+ 2d'^ + d) + 4r2 = 0. (VII.35) 

Upon eliminating c and d from these equations, we get 

j,‘2 £ff £t o 

- - - - -4/r(c' + 2d') — -r‘^f{c' + 2d') — Afr^c'd' + 8r^ + 2/^(d)^ 

- r^f{c'^ + 4d'2) - fr\c" + 2d") + Al^cd - 2/ = l^Wu (VII.36) 


rf" + 3/' + 3r/'c' + 2rf'd' + 2fr{df + 8/d' + 4r/(d')2 + 2/rc" + 4/rd" 
-8r = -2rfWrr (VII.37) 

At the second order flVII.361) & flVII.371) take the form 


- 4/r(c'2 + 2d'f) - ^r^f'ic', + 2d',) - fr\c'' + 2d'') = (VII.38) 

rV'(3c'2 + 2d',) + 8rfd', + 2fr\c'' + 2d'') = -2«:rV(<^>'i)' (VII.39) 

Let C 2 + 2d2 = X. Then flVII.38p becomes 


-AfrX' - -r^fX' - fr^X" = kI\^i)^ 


x„ ^ !+±ir+A" = -k;- 


(* 1 )= 


f^2 

X" + P{r)X' = -kI 


fr^ 


/^2 


where P{r) = 


ifr+^Pf 
fP • 


Therefore, 




X' = -nde-rnr)<ir 


(VII.40) 
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Thus 


c 


/ 

2 


1 

Arf — 2f'r‘^ 




/"(<i>i)2)-X'(4r/ + 2rY) 


(VII.41) 


4 = (VII.42) 

By integrating both sides of (1VIL41I1 and (1VII.42I1 we can obtain C2 and d 2 ■ 
Putting 02 = &2 = 0 in (111.181) we get the equation governing the dynamics 
of the system at the third order, which is : 


+ L<h3 = ^(C2, da, <hi) (VII.43) 

where S = — (c’a + 2^2)*hi + ^{c '2 + 2^2) 


VIII Appendix B 

The equations of motion for the Einstein-Gauss-Bonnet (EGB) scalar field 
system are: 

Gab + J^Qab + hHab = SllGidaCpdbCt) - ^gab{d(j)Y) 

9^^^aSb<t> = 0 (VIIL44) 

where 

Hab = 2{RRab - 2Rafj,R\ - 2R^'"R^aub + Ra'^'^Rbiiup) 

~gab{R^ - ^R^'^Rpu + R^'^^^Rpuap) (VIIL45) 

A3 is the Gauss-Bonnet coupling. As shown in |22], this can be cast in hrst- 
order form as in the Einstein gravity-scalar held case. Studying the system 
order by order in weakly nonlinear perturbation theory as in II. 1, there are 
corrections to the corresponding equations in Einstein gravity beyond the 
linear order in e. The equations to O(e^) are: 

At the hrst order: (the equation is identical to the Einstein gravity case) 

- 01 + \ da^itSiV? x(0a:0l)) = 0 (VIII.46) 

tan'^ X 
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At second order: 


^2 = _^2A3) Jo (VIIL47) 

= (, + {4>if)dy (VIII.48) 

(1 - 2A3) sin (x) Jo 

At the third order: 

~ 03 + 7— 3 —'9a;(tan^ x{dx4>3)) = 20i((52 + A2) + 4)1(62 + A2) + 01 (A2 + <^2) 
tan'* X 

(VIII.49) 

At the fourth order: 


A4 = 


COS^ X 
sin^ X 



—^3 

-(I- 2 A 3 ) 


(^)’ + (r=W ^ 

+ 252 ( 0 i )^ + ^2(01)^ — A 2 ( 0 i)^) dy^ 

(VIIL50) 


(1 — 2A ^ (20103 + 20103 + 2^2(01)^ + 2^2(0 

2A3 (^2^2 

(1 - 2A3) sin^iy) 


dy (VIIL51) 


Let us now consider the following rescaling of the held: 


0 ^ , 

(1-2A3)V2 ^ 


(VIII.52) 


At the hrst order, there is no difference in the equation for 0 and 0. 

- 01 + \ 0a;(tan^ x(dxipi)) = 0 (VIII.53) 

tan"’ X 

At the second order, this equation is identical to the Einstein gravity-scalar 
held system for 0. 

62 = - [ cosi/sini/((00^(0i)^)(ii/ (VIII.54) 

Jo 
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(VIIL55) 


sin^ X Jq 

At the third order; 

~ t/’s + 7-3— dx{t&n^ x{dx1p3)) = 2'^i(<52 + ^ 2 ) + t/^l(52 + A 2 ) + t/'l(^2 + ^ 2 ) 

tan'^ X 

(VIII.56) 

At the fourth order, there are terms proportional to A 3 that cannot be scaled 
away. 


A4 = 


COS^ X 
sin^ X 



- -A 3 / (Aa)^ y 

-(1 - 2 A 3 ) Vcos^(|/)/ 


+ tan^ yi^ip'iip's + 2'ipi'ip^ + 252{'ipif 


+A2{:ipif 


dy} 


(VIII.57) 


<54 = 



cosy siny (2^7 '^/>3 + 2 ^/>iV ^3 + 2 A 2 (V^i)^ + 252{'<piY) 


2A3 S2A2 
(1 - 2 A 3 ) sin2(y) 


jdy (VIII.58) 
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